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We address the question whether features known from quantum chromodynamics (QCD) can possibly also 

■ - - ■ show up in solid-state physics. It is shown that spinless fermions of charge e on a checkerboard lattice with 

^M ' nearest-neighbor repulsion provide for a simple model of confined fractional charges. After defining a proper 

^^ I vacuum the system supports excitations with charges ±e/2 attached to the ends of strings. There is a constant 

^^ , confining force acting between the fractional charges. It results from a reduction of vacuum fluctuations and a 

' polarization of the vacuum in the vicinity of the connecting strings. 

^ ■ 

•^ , One of the most fascinating aspects of quantum chromodynamics (QCD) is the existence of particles with 

fractional charges, i.e., of quarks and their confinement''^. Of particular interest is that the confining force is 

constant, i.e., independent of the distance of the separated particles. The question may be posed whether this 

physical feature is restricted to QCD or whether it occurs in a similar form in other areas of physics, e.g., in solid- 

~— 1 state theory. There are well known examples where parts of the above mentioned physics show up in solids. We 

q_) mention here that in heavily doped polyacetylen excitations with fractional charges may exist^. In two dimensions 

J the fractional quantum Hall effect is another example for the appearance of fractionally charged excitations'*. But 

in both cases there is no constant confining force. On the other hand, confinement is nearly realized, e.g., in an 

Ising antiferromagnet with a hole put into it. When this hole is moving with hopping matrix element t it generates 

"j^ ' a string of frustrated spins and therefore a constant confining force - almost, one should add, because there exist 

Ch so-called Trugman paths^ which enable the hole to avoid generating such string. The probability for those paths 

I is small, though, i.e., the effective hopping matrix element is of order t^ /J'^ where J is the coupling constant of 

■ O interacting neighboring Ising spins. But in this example there are no fractional charges appearing. 

5 ' The aim of the present communication is to point out a simple Hamiltonian for spinless (or fully spin polarized) 

fermions of charge e defined on a frustrated lattice, here a checkerboard lattice, which leads to excitations with 

fractional charges ±e/2 and a constant confining force. The advantage of studying such a model system is that it 

enables us to describe in detail the origin of the confining energy and constant confining force. We consider the 

^"^ following Hamiltonian on a checkerboard lattice^ 

^' H = -tY,(c]cj+h.c\+vY,ninj (1) 

^ , ('7> (ij) 

^' t t 

\^ \ where c\ {ci) create (destroy) spinless fermions on site / and n, = c\ci. The notation {ij) refers to neighboring sites. 

^D . The checkerboard lattice can be considered as a projection of a pyrochlore lattice onto a plane. The latter consists 

of corner sharing tetrahedra and therefore each of its sites has six nearest neighbors. We are interested here in the 

limit y> |f|. 

First we want to define a vacuum state for our system. For that purpose we consider a lattice of A^ sites with N /2 

, , particles. When f = the repulsion energy is minimized if every crisscrossed square is occupied by two particles 

r^ ("tetrahedron" rule)^. In that case the (classical) ground state is approximately (5)^^ - fold degenerate as known 

^ from the ice problem^. One special of those configurations \^\ ) is shown in Fig. la. Here we have connected 
. . neighboring occupied sites by a solid line. Assuming periodic boundary conditions those lines form loops which 
are either contractible or go around the torus. The degeneracy of the ground states is partially lifted by the hopping 
term^. To lowest order in t/V this occurs through ring hopping as described by the effective Hamiltonian 
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^eff = -r7rL444^^5^^3Cji • (2) 

It acts on the subspace spanned by the different ground-state configurations and within that space on hexagons 
with sites ji ,...J6 which are alternating occupied and empty. One of those hexagons is indicated in Fig. la. Note 
that ring exchange is quite common in He^ (see Ref.''') and in frustrated spin systems". When matrix elements 
((|)J |//eff|(|) ■ ) are calculated it matters for the sign of the result whether the site in the center of the hexagon is 
occupied or empty. For that reason it is advantageous to rewrite //gff in the form 
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Figure 1: (a) Ground-state configuration maximizing the number of flipable hexagons, one of which is explicitely shown. 
Neighboring occupied sites (dots) are connected by a solid line, (b) Two defects with charge e/2 of an arbitrarily chosen 
allowed configuration are connected by a string consisting of an odd number of particles, shown by the black line. 



//,., 



■L 



<X>)(<X>I l<X>)«X>l 



H.c 



(3) 



with g = 12t^/V^. The sum is over all horizontal and vertical hexagons of the system. The remaining notation 
is self-explanatory. Note that terms of order t^/V yield identical contributions for all ground states and therefore 
do not lift their degeneracy. Exact diagonalization of finite clusters with different geometries shows that within 
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the quantum-mechanical ground state |v|/o) = LiCCi|<l'; ) those |(|)J ) will have largest amplitudes a, which contain 

the largest numbers of flipable hexagons. An investigation of the configurations |(|)| ) of finite clusters with up 
to 1000 sites shows that the ones with the largest number of flipable hexagons form squiggles. The latter have 
been first identified by Penc and Shannon'^. Their unit cell has 40 sites. From translation and rotation of the 

unit cell a corresponding number of different configurations |(|)J ) can be generated which maximize the number 
Nfi of flipable hexagons. One of them is shown explicitely in Fig. la. These configurations correspond to two 
different charge orderings which are related by a rotation of 90° around a symmetry axis. This is seen by counting 
the number of particles along the lines denoted by 1 and 2 in Fig. 3a. The charge ratio in the two directions is 
{n^) : («y) =2:3 and for the rotated configuration 3 : 2. 

When in a small cluster with periodic boundary conditions the squiggles can not develop, configurations with 
maximum Nf[ have a different shape. For quantum mechanical calculations we have used a cluster of size \/72 x 
\/72 with periodic boundary conditions. From the numerical analysis of such a cluster we have determined the 
distribution of the weight of configurations with given Nf^. The result is shown in Fig. 2. A single configurations 

\(^j ) contributes to |\|/o) the more, the more flipable hexagons it has. But due to their large number, the largest 

total contribution to the ground state of, e.g., a 72 site cluster comes from configurations |(|)- ) with 14 flipable 
hexagons, i.e., N^ = 14. Since the squiggle phase cannot be realized in the considered cluster, we find a ratio 
(%) : (rty) = 1:2 and 2 : 1 instead for the configurations which maximize N^. The Hamiltonian (3) conserves that 
ratio. Since charge ordering is solely caused by the effective kinetic Hamiltonian (3), we have here a nice example 
of order by disorder. The quantum mechanical ground state is therefore 2-fold degenerate. We define the vacuum 
by choosing one of the ground states, i.e., we are dealing with a symmetry broken (degenerate) vacuum. 

Next we place an additional particle onto an empty site of the checkerboard lattice. Hereby we select that 
particular ground state which in the presence of the added particle gives the lowest energy. As a consequence two 
neighboring tetrahedra have three particles attached to them. When the hopping term of Eq. (1) is applied the two 
tetrahedra separate. That does not cost any additional repulsive energy. Each of the two tetrahedra has a charge e/2 
attached to it, if the charge of the added particle was e. Separating the two tetrahedra with three particles generates 
a string consisting of an odd number of sites which is linked to two loops with an even number of sites each. This 
is shown in Fig. lb for an arbitrarily chosen configuration. The fractional charges sit at those linking points. Next 
we determine the changes of the kinetic energy density due to the presence of the two charges e/2. This is done by 
keeping the two tetrahedra with the fractional charges fixed at and r and determining the ground state \(/o(0,r) 
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Figure 2: Average weights of configurations with Nf[ flipable hexagons in the quantum mechanical ground state of a 72 site 
checkerboard cluster. 



and its energy for a cluster of 72 sites. The latter can be decomposed into local contributions by calculating the 
following expectation value of the energy at site ; 

B,- = -f E ('^o(0,r)|(|<^)(<^|-|<^)«^|+H.c.)|^o(0,r)) . (4) 

The sum is over all hexagons containing site /. The result is shown in Fig. 3a,b. One notices a decrease in 
kinetic energy in the region between the two fractional charges, i.e., along the connecting string. The reason is 
that due to the topological changes caused by the generated string the number of flipable hexagons positioned 
between the fractional charges is decreased. This energy decrease is proportional to the length of the generated 
string and implies a constant confining force. The origin of that force is clear. It is due to a reduction of vacuum 
fluctuations caused by the topological changes in the vacuum when the fractional charges are separated. Consider 
a configuration with two defects of charge e/2 separated along a diagonal (see Fig. 3b). Then we find that the 
maximal number of flipable hexagons is reduced by 2d, where d is the number of crisscrossed squares between the 
defects. We have also determined the density changes caused by the separation of the fractional charges (see Fig. 
3c,d). They are obtained from the ground state in the presence of the two defects |\(/o(0,r)) by simply calculating 



5«,- = (vi/o(0,r) I Hi I \i/o(0,r)) - (\(/o | «,■ | \|/o) 



(5) 



While the total changes add up to zero, the vacuum is modified by the breaking up of the charge. It is polarized 
along the connecting string. The above considerations can be also formulated in terms of a dimer model'^. For 
that purpose one connects the centers of the crisscrossed squares ending up with a square lattice. If the site on a 
link is occupied a dimer is attached to it. The tetrahedron rule implies that each site of the square lattice is touched 
by two dimers. 

The same findings discussed here apply also when a particle is removed from the ground state or when a particle- 
hole excitation is generated out of the ground state. They remain nearly unmodified for the case of quarter filling. 
We also expect to find a similar picture for the kagome lattice, another frustrated lattice in which the classical 
ground-state manifold can be described by a dimer model on a bipartite lattice. We have not discussed here the 
three dimensional pyrochlore structure where also fractional charges may occur^ because the numerics is much 
more involved. Neither have we discussed possible experimental verification of fractional charges in frustrated 
lattices^^, e.g., in the rather common spinels (for the latter see, e.g., Ref.^^). Instead we refer to the literature. 

Summarizing, we have found that a simple Hamiltonian ( 1 ) when applied to a half-filled frustrated lattice leads in 
the strong coupling limit to fractional charges ±e/2 with a constant confining force between them. The latter is due 
to modifications of the symmetry broken vacuum between the charges, i.e., the vacuum fluctuations are diminished 
and the vacuum is polarized along the string connecting the fractional charges. Our findings suggest that a number 
of features known from QCD are also expected to occur in solid-state physics. Conversely, one would hope that by 
studying frustrated lattices or dimer models one might be able to obtain better insight into certain aspects of QCD. 
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Figure 3: (a) and (b): local loss of kinetic energy due to a separation of two fractionally charged defects. Radii of the dots are 
proportional to the local energy loss. The positions of defects are shown by red colored crisscrossed square, (c) and (d): red 
and blue dots show an increase (decrease) of the local density (vacuum polarization). 

A crucial difference between the confinement in the introduced model and in QCD is the fact that we can adjust 
the confining force by varying the ratio tjY. For the considered limit |f | ^ V, the fractionally charged particles 
form bound pairs of large spatial extent'^. Thus we expect that charge fractionalization has observable effects even 
if the fractionally charged particles are confined. For instance an exchange of fractionally charged particles might 
occur at finite doping. 
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